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Abstract 

Using the published Kamiokande data of the multi-GeV atmospheric neu- 
trinos, we have searched the optimum set of the neutrino oscillation parame- 
ters among three flavors. It is found that \ 2 is minimized for (Am^, Am^) 
= (3.8xl(T 2 eV 2 , 1.4xl0- 2 eV 2 ), (6 l2 , 6 13 , 6 23 ) = (19°,43°,41°) with X min 
= 3.2 (42%CL). The sets of parameters (Am 2 21 , Am 2 x ) = (e>(KT n eV 2 ) or 
O(10 _5 eV 2 ), (9(10 _2 eV 2 )) which are suggested by the two flavor analysis fall 
within 0.7cr. 
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There has been much interest in atmospheric neutrinos |!J 01 @ @] > 
which might give us an evidence for neutrino oscillations. While NUSEX [|J 
and Frejus || have reported consistency between the data and the predictions 
on atmospheric neutrino flux || [0], Kamiokande, 1MB and Soudan-2 have 
reported discrepancy. In particular, the Kamiokande group claimed that 
their multi-GeV data suggests the mass squared difference of neutrino is of 
order 10" 2 eV 2 [g. 

People have studied neutrino oscillations among three flavors ||, and it 
has been shown recently |10| that the mass squared differences and the mixing 
angles have strong constraints from various experiments. The analysis of 
the multi-GeV atmospheric neutrino data by the Kamiokande group was 
based on the framework of neutrino oscillation between two flavors^], and it 
is important to see what happens if we analyze the data in the three flavor 
framework. In this paper we will analyze the published multi-GeV data [g 
of the Kamiokande atmospheric neutrino experiment, taking into account 
mixings among three flavor neutrinos. Unlike other works in Ref. [[RJ, we 



will take the matter effect [13] into consideration, and evaluate the number 

of events by summing over the energy and the zenith angle of neutrinos, to 

reproduce the original analysis by the Kamiokande group as much as possible. 

Throughout this paper we will restrict our discussions only to the multi-GeV 

data by the Kamiokande group, not only because the Monte Carlo result for 

the neutrino energy spectrum is available only in Ref. 0, but also because 

this is the only data which gives both the upper and the lower bound on the 

mass squared difference of neutrinos. 

We start with the Dirac equation for three flavors of neutrinos with mass 

2 A couple of works have discussed Kamiokande's analysis from the viewpoints 

which are different from ours. 
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in matter [fT5 |: 
• d T / n 



C/diag (Ei, £ 2 , £ 3 ) + diag (A(x), 0, 0)1 *( 



(1) 



Here i?j = yp 2 + m| is the energy of the neutrino, \fr(x) = (v e (x), u^x), v T {x)) 7 
is the wave function of the neutrinos in the flavor basis, A(x) = \f2GpN e (x) stands 
for the effect due to the charged current interactions between v e and electrons 
in matter [fT5 |. 



U 




C12C13 
-S12C23 - C12S23S13 

\ S12S23 — Ci 2 C 2 3Si3 



S12C13 
C12C23 - Si 2 S 2 3Sl3 
— C12S23 — S12C23S13 




(2) 



with c. 



1.1 



cos 6*. 



sin fly is the orthogonal mixing matrix of neutrinos, 
and we will not discuss the CP violating phase of the mixing matrix here for 
simplicity []. 

The number of the expected charged leptons £ a ( £ a = or r ) with 
energy q from a scattering v a e — > v e i a is given by 



dE / dQ 



<7max 



= n T £ 

(3=e,fi 



u a ;E,Q) (a = e,/i) (3) 



Here Fp(E,0) is the flux of atmospheric neutrino z/g with energy from 
the zenith angle G, ut is the effective number of target nucleons, e(q) is 
the detection efficiency function for charged leptons £ a , da a (E,q)/dq is the 



differential cross section of the interaction v n e 



v, 



[a 



eor /i), P(u 



3 Even if we include the CP violating phase 8 of the mixing matrix, the effect of S always 
appears in the combination of si^e s . S13 has to be small because of the constraints from 
the reactor experiments and from the solar neutrino observations, as we will discuss below. 
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is a ;E,Q) = 1(^(0) |z/ Q (L))| 2 is the probability of Vp — > v a transitions with 
energy E after traveling a distance 



(R + h) 2 - R 2 sin 2 G-R cos 6, (4) 

where R is the radius of the Earth, h ~15Km is the altitude at which atmo- 
spheric neutrinos are produced. 

To reproduce the analysis of the multi-GeV data by the Kamiokande 
group, one needs the quantity 

fpa(E, 9) = n T jf*"" dq e(q)F p (E, 6) g) (a,(3 = e,fi) (5) 

for each E and O, which is not given in ||. However, the quantity 

gp a (E)= r dB Ua(E,e), (6) 
Jo 

which is obtained by integrating ([5]) over B, is given in the Fig.2 (d)-(f) in 
Ref. 0. The zenith angle dependence np{E, 0) of the atmospheric neutrino 
flux for various neutrino energy E has been given in Ref. |7J in detail. Here 
we multiply the quantity gp a (E) by the zenith angle dependence in Ref. [0 
with suitable normalization, and adopt the quantity 

7 I jp n\ _ 9Pa(E)rip(E, 6) ( o \ (n\ 

Ua (E,Q) = — M ^ Q) (a,(3 = e,ri (7) 

instead of the original quantity fp a (E,Q) used in 0. (0) is the important 
assumption of the present analysis. (|^) is not exactly the same as fp a (E,Q) 
in the original analysis @, but this is almost the best which can be done 
with the published data in [§]. 

We have solved (p} numerically for each E (1(T 1/20 GeV < E < 10 2 
GeV) and evaluated the number of events for a given range of the zenith 
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angle = cos" 1 ^) < < J+1 = cos"^^) (1 < j < 5) 



Xj 4 = N(fi, Qj<Q< 

= (l + a )(l + 9 / d6 / dE 

f^(E, e)P(u tl - »/ M ; E, 6) + / eM (£, 0)P(z/ e - ^; £?, 9) 

xj = N( e , e, < e < e i+ i) 

= (l + a )(l-£) / de / dE 

ffie{E, e)P(u lx - i/ e ; -E, 0) + f ee (E, Q)P(u e -> z/ e ; £, 0)] . (8) 

Several groups || have given predictions on the flux of atmospheric neutri- 
nos but they differ from one another in the magnitudes, and the Kamiokande 
group assumed that the errors of the overall normalization 1 + a and the rel- 
ative normalization 1 + 0/2 are a a =30% and a p=12%, respectively. Here we 
regard these factors a and as free parameters of the theory, and adopt the 



following x |M 



a=e,fj,j=l \ 



Nf In -i- 

, at;, 



(9) 



where A 7 "" (a = e, /i; 1 < j < 5) is the data for each zenith angle 0j < < 



Oj+i. The theoretical prediction X? 



e, (/,;!< j < 5) depends on seven 



free parameters (Am^Am^; ^12^13,^23; ck,/?) , where Am^ = mj — m|, so 
(Q) is expected to obey a x 2 distribution with 10—7=3 degrees of freedom. 
The number of degrees of freedom in the present analysis is smaller than the 
original one by the Kamiokande group (5x8+5—2=83). 

We could not reproduce exactly the zenith angle distributions in Fig. 3 
in Ref. ||. In particular our prediction for e-like events near cos© ~ —1 has 
a larger difference with the data than Kamiokande's does, and this difference 
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seems to be important to discuss the magnitude of x 2 laterQ. Presumably this 
discrepancy arises not only because the data that we are using is different 
from the original one in Ref. 0, but also because the Kamiokande group 
has taken into account the smearing effect on the resolution of the angle 



(15° ~ 20°) and the effects of backgrounds [15|. Throughout this paper 
we discuss the goodness of fit and the confidence level of set of the parameters 
etc. based on our calculation with (0). 

The value of x 2 is affected to some extent by the presence of matter, 
and it is necessary to take into consideration the contribution of the second 
term in ([!]). Evaluation of x 2 requires a lot of CPU time of a computer since 
one has to solve ([I]) numerically for each E and and plug it into @. We 
have meshed each parameter region into ten points (Am?- = 10~ 5+ ^ 2 , % = 
£7r/20 (0 < £ < 10)) and the evaluated the value of x 2 ■ Furthermore, using 
the gradient- search method described in Ref. fll6| , we have found that x 2 
has the minimum value for 

(Am^Am^) = (3.8 x 10~ 2 eV 2 , 1.4 x 10~ 2 eV 2 ) 
(012,013,023) = (19°, 43°, 41°) 

(a, /3) = (2.8 x 10~\ -5.0 x 10~ 2 ) (10) 

with Xmm = 3-2. Note that the deviation of the two normalization factors 

1 + a and 1 + (5/2 from unity is within the errors Q =3O% 073=12% assumed 

in Ref. 0. We have also calculated the value of the modified chi square 

X 2 = x 2 + 0/2 l a a + f3 2 / a 1 with the weight for the errors a a =30% 073= 12%, 

and we have found that the conclusions in the following discussions do not 

change very much with x 2 instead of x 2 - 

4 If we try to fit the data with only two parameters a and (3 as in Ref. p2f , then the 
minimum value of x 2 is 21, which suggests that the 9 independent solution is excluded at 
the 99% confidence level in our analysis. 
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The zenith angle distributions of the e-like events, the /i-like events and 
the double ratio R = (/V e )data/(AV e )MC are §i ven i n Fig.l an d Fig.2 for the 
optimum set of parameters. 



(Insert Fig.l and Fig. 2 here.) 

The degrees of freedom of our analysis is 3, so the value of the reduced chi 
square is 1.1, which corresponds to 42 % confidence level. This suggests that 
our fit in the present analysis is not particularly good, but as we mentioned 
earlier, this is probably due to the fact that the data from which we start is 
poorer than the original one by the Kamiokande group 

The region of the parameter space which is allowed at 90 % CL is given 
by x 2 — Xmin + 12 for seven free parameters. However, because it requires 
a lot of CPU time of a computer to solve ([!]) numerically, we could not give 
sets of contours in the parameter space. In fact, since the value % ~ Xmin is 
rather large, if we project the allowed region onto the Am^-Am^ plane, it is 
conceivable that we have disjointed regions on this plane, and the calculation 
would be extremely tedious. So we restrict our analysis to a special case of 
particular interest here. 

We have evaluated x 2 f°r the sets of parameters, which are suggested by 
the solutions for the solar neutrino problem [|13j [[H]]. In order not to spoil 
the success of these scenarios |13| (T7j based on the two flavor framework, 
we consider only the case in which \U e s\ is small, since we have the formula 
which relates the probability P^\v e — > u e ) in the three flavor analysis to 
P {ye — ¥ v e) in the two flavor one [|T8[|: 

P (3) (z7e ^ Ve . A{x)) = (1 _ \U*\*)*pV\v e - u e ; (1 - \U e3 \ 2 )A{x)) + \U e3 \\{ll) 

We note in passing that the constraint for \U e s\ also comes from the reactor 
experiments [BJ|, which suggests \U e3 \ 2 = s± 3 < 10 _1 or 1 — \ U e 3\ 2 = c 2 3 < 10 _1 
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for AmJi ~ a few lCT 2 eV 2 . 

Irrespective of whether we consider the vacuum solution (Am^ ~ (9(10 _11 eV 2 )) 
or the MSW solution (Am^ ~ (9(10~ 5 eV 2 )) for the solar neutrino, the mass 
squared difference AmJi is negligible compared to the contribution of the 
matter effect A(x) in (j^) and the other mass squared difference Am^, which 
should be at least of order 10 _2 eV 2 to account for the zenith angle depen- 
dence of the Kamiokande multi-GeV data. Thus we consider the case 

(Am2 1 ,sin 2 2^ 12 ) = (Am 2 , sin 2 26) Q 

' (C(10- n eV 2 ),O(l)), (vacuum solution) 

= < (C(10- 5 eV 2 ), 0(l(r 2 )), (small angle MSW solution) 

^ (O(10- 5 eV 2 ), 0(1)), (large angle MSW solution) 
< 13 < 5° 

Anijj , 622, = arbitrary (12) 
With this constraint we have found that x 2 has the minimum value for 

Am 2 31 = 3 x l(T 2 eV 2 

9i3 = 5° 

9 23 = 40° or 50° 
(a, (3) = (3.6 x 10 _1 ,-3.7 x 1(T 2 ) (13) 

with 

X 2 = 9.6. (14) 

Deviation of each parameter in this case is (x 2 ~ Xmm)/7 = 0-9> so we conclude 
that this set of parameters falls within 0.7a for all three cases in ([12]). In 
fact we observe that any set of the parameters with the constraints falls 
within 0.7a as long as Am^ <C Am 2 2 < Am^ is satisfied (arbitrary # 12 is 
allowed in this case). In ( |T3"D the error of a is a little too large compared to 
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what the Kamiokande group assumed, but even if we take a = 3.0 x 10 _1 
with all other parameters the same as in (|T3|), we find that the solution falls 
within 0.8a. The reason that we have weaker constraints in this analysis 
than in Ref. || is because we have larger numbers of free parameters, but 
this is inevitable as long as one assumes the general mixings among three 
flavors of neutrinos. 

In this paper we have analyzed the multi-GeV atmospheric neutrino 
data by the Kamiokande group based on the framework of three flavor 
neutrino oscillations, and have shown that the best fit is obtained for the 
set of parameters Am^ ~ Am 2 ^ ~ C (10 _2 eV 2 ). We have also shown 
that the popular set of parameters ( (Am^, sin 2 26*12) = (Am 2 , sin 2 29) , 
(Am 2 1 ,sin 2 2^ 13 ) = (£>(10~ 2 eV 2 ), 0(1)), 13 ~ ) fall within 0.7a. The 
minimum value of x 2 is 3.2 for 3 degrees of freedom, and the fit based on the 
hypothesis of neutrino oscillations is not particularly good due to that fact 
that we used only the information published in Ref. 0. We hope that the 
situation will be improved much more when the SuperKamiokande experi- 
ment starts. If we combine the results here with other experimental data, 
then we get even stronger constraints, which will be reported somewhere jOD] . 
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Figures 



Fig.l (a),(b) Zenith angle distributions for the e-like and /z-like multi-GeV 
events. The squares with error bars are data and the histograms stand 
for the predictions without neutrino oscillations (solid lines), and with 
neutrino oscillations (dashed lines for the optimum set of parameters 
(AmL AmiO = (3.8xl(T 2 eV 2 , 1.4xl(r 2 eV 2 ), (0 12 , 9 13 , 23 ) = 
(19°, 43°, 41°)), respectively. These quantities are obtained by multi- 
plying the values in Fig. 3(d) in Ref. @ by the zenith angle dependence 
of the flux in Ref. 0. 

Fig. 2 Zenith angle distribution of the double ratio R = (/V e )data/(/V e )MC- 
The solid lines stand for the prediction with neutrino oscillations for 
the optimum set of parameters. All the quantities are calculated based 
on the same assumption as in Fig.l. 
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Fig.1(a) 
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